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1 Introduction 



Numerous observations of SN 1987 A suggest that extensive mixing has taken place in the exploding 
envelope of the progenitor star Sk -69 202. The early detection of X and 7-rays, the broad profiles 
of infrared Fe and Co lines, as well as the shape of the light curve cannot be explained without 
assuming that clumps of newly synthesized 56 Ni, from layers close to the collapsed core, have 



penetrated into the hydrogen envelope (see the reviews of Arnett et al. 198E and Muller 1998, 
and the references therein). That such mixing is probably generic in core collapse supernovae 
is indicated by spectroscopic studies of SN 1987 F, SN 1988 A, SN 1993 J ( |5pyromilio, 1994 and 
references therein) and SN 1995 V (Fassia et al., 1998). Furthermore, it might also explain the 



detection of fast moving clumps of metal-enriched material in the Vela (Aschenbach et al., 1995), 
Cas A (Anderson et al., 1994) and Puppis A ( Winkler fe Kirshncr, 1985| ) supernova remnants 
as well as the isotopic composition of specific SiC grains with possible supernova origin found in 
primitive meteorites (Nittlcr et al. 1996; Amari this volume). But even more important from the 
point of view of supernova modellers is the fact, that a detailed understanding of the problem 
of nucleosynthesis and mixing can give us invaluable information about the explosion mechanism 
itself. Thus, the observations have instigated theoretical work on multidimensional supernova 
models which focused eithe r on the role of convection occuring within the fir s t second of a delayed , 



1998; Janka & Muller, 1996; Burrows et al., 199E 



Hcrant et al., 1994; 


Miller et al., 1993 


), or on the growth of Ray lcigh- Taylor instabilitites during 


the late evolutionary stages ( 


Nagataki et al., 1998: 


Hcrant & Benz, 1992; 


Muller et al., 1991; 



Yamada fe Sato, 1991] ; [fiachisu et al., 1990| ). 

However, multidimensional simulations which follow the evolution of the stalled supernova 
shock from its revival due to neutrino heating, until its emergence from the stellar surface have 
not yet been performed. Due to the presence of vastly different spatial and temporal scales 
and the range of physical processes involved during the early stages following core collapse, all 
studies of mixing in core collapse supernovae have hitherto neglected the influence of neutrino- 
driven convection in seeding the Rayleigh- Taylor instabilities. Instead, a shock wave was created 
artificially by depositing the explosion energy near the center of a pre-collapse progenitor model 
and following the propagation of the shock in one spatial dimension until it had reached one of 
the unstable composition interfaces. This was either chosen to be the He/H interface in case of 
Type II supernova models (e.g. Hcrant & Benz 1992; Muller et al. 1991) or the C+O/He interface 
in case of Type lb models (Hachisu et al., 1994). Only then were the ID models mapped to a 2D 
grid and the rest of the evolution followed with a multidimensional code. A somewhat different, 



two-dimensional approach, was chosen in the recent calculations of Nagataki et al. (1997) who 
initiated the explosion using a parameterized, aspherical shock wave and computed the resulting 
nucleosy nthesis using a mark er particle approximation. These models have been subsequently 
used by Nagataki et al. (1998j ) for a study of Rayleigh- Taylor instabilities at the He/H interface. 
Still, however, these calculations do not address the complications introduced by the explosion 
mechanism and thus suffer from a number of assumptions. Furthermore, their numerical resolution 
appears to be hardly sufficient to resolve instabilities which occur within the first minutes of the 
evolution. 



2 



In the present contribution a first step towards a more consistent multidimensional picture of 
core collapse supernovae is attempted by trying to answer the questions 



What happens in the first minutes and hours of a core collapse supernova if neutrino heating 
indeed succeeds in reviving the supernova shock? 



• Can neutrino-driven convection in conjunction with the later Rayleigh- Taylor instabilities 
lead to the high iron velocities observed in SN 1987 A? 

For this purpose we have carried out high-resolution 2D supernova simulations which for the 
first time cover the neutrino-driven initiation of the explosion, the accompanying convection and 
nucleosynthesis as well as the Rayleigh- Taylor mixing. In the following, we present preliminary 
results from these calculations, focusing on the first ~ 300 seconds of evolution. A summary of 
our work can also be found in Kifonidis ct al. (1999] ). 



2 Numerical Method and Initial Data 

We split our simulation into two stages. The early evolution (t < Is) which encompasses shock 
revival by neutrino heating, neutrino-driven convection and explosive nucleosynthesis is followed 
with a version of the HERAKLES code (T. Plewa & E. Miiller, in preparation). This hydrody- 
namics code solves the multidimensional hydrodynamic equations using the direct Eulerian version 



of the Piecewise Parabolic Method (Colella & Woodward, 1984) augmented by the Consistent Mul- 



tifluid Advection (CMA) scheme of Plewa & Miiller (1999) in order to guarantee exact conservation 



of nuclear species. We have added the input physics (neutrino source terms, equation of state, 



boundary conditions, gravitational solver) described in lanka & Miiller (1996) (henceforth JM96) 



with the following modifications. General rclativistic corrections are made to the gravitational po- 



tential following Van Riper (1979). A 14-isotope network is incorporated in order to compute the 
explosive nucleosynthesis. It includes the 13 a-nuclei from 4 He to 56 Ni and a representative tracer 
nucleus which is used to monitor the distribution of the neutrino-heated, neutron-rich material 
and to replace the 56 Ni production when Y e drops below ~ 0.49 (cf. Thielemann et al., 1996 ). 



We start our calculations 20 ms after core bounce from a model of Bruenn (1993) wh o has 
followed core-collapse and bounce in the 15 M Q progenitor model of |Woosley et al. (1988 ). The 



model is mapped to a 2D grid consisting of 400 radial zones (3.17 x 10 6 cm < r < 1.7 x 10 9 cm), 
and 180 angular zones (0 < 9 < ir; cf. JM96 for details). A random initial seed perturbation is 
added to the velocity field with a modulus of 10 -3 of the (radial) velocity of Bruenn's post-collapse 
model. The calculations are carried up to 885 ms. At this time the explosion energy has saturated 
and essentially all nuclear reactions have frozen out. We will henceforth refer to this calculation 
as our "explosion model" . 

The subsequent shock propagation through the stellar envelope and the growth of Rayleigh- 
Taylor instabilities is followed with the AMRA Adaptive Mesh Refinement (AMR) code (T. Plewa 
& E. Miiller, in preparation). This code uses a different variant of HERAKLES as its hydrody- 
namics solver which does not include the neutrino physics. Gravity is neglected in the AMR 
calculations since, as is the case for the neutrino source terms, it does not influence the propaga- 
tion of the shock during late evolutionary stages. However, gravity is important for determining 
the amount of fallback, a problem which is outside the scope of the present study. The equation of 
state takes into account contributions from photons, non-degenerate electrons, e + e _ -pairs, 1 H, and 
the nuclei included in the reaction network. The AMR calculations are started with the inner and 
outer boundaries located at r- m = 10 8 cm (i.e. inside the hot bubble containing the neutrino-driven 
wind) and r out = 2 x 10 10 cm, respectively. No further seed perturbations are added. We use up 
to four levels of mesh refinement and refinement factors of four in each grid direction, yielding a 
maximum resolution equivalent to that of a uniform grid of 3072 x 768 zones. In order to keep 
the radial resolution as high as possible during any given evolutionary time, we do not include the 
entire star but allow the code to expand the radial extent of the base grid by a factor of 2 to 4 
whenever the supernova shock is approaching the outer grid boundary. The latter is moved from 
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Figure 1: Evolution of the entropy in our explosion model within the first 300 ms after the start 
of the calculations. Neutrino-heated, deleptonized material rises in large bubbles and distorts the 
shock (outermost discontinuity). 

its initial value out to r out = 1.1 x 10 12 cm at t = 300s. Reflecting boundary conditions are used 
at 9 = and 8 = tt and free outflow is allowed across the inner and outer radial boundaries. 



3 Nucleosynthesis and Neutrino-Driven Convection 

The general features of our explosion model are comparable to the models of JM96. For the 
initial neutrino luminosities, which are prescribed at the inner boundary, somewhat below the 
neutrino sphere, and decay with time as described in JM96, we have adopted a value of lP v = 
2.8125 x 10 52 erg/s for the electron neutrinos and lP v — 2.375 x 10 52 erg/s for the heavy lepton 



neutrinos (with v x 



,v^,v T ,v T ). The parameters describing lepton and energy loss of the inner 



iron core were set to AY; = 0.0875, Ae = 0.0625 (cf. JM96). The neutrino spectra are the same 
as in JM96. 

For the chosen neutrino luminosities the shock starts to move out of the iron core almost 
immediately. Convection between shock and gain radius sets in ~ 30 ms after the start of the 
simulation in form of rising blobs of neutrino-heated, deleptonized material (with Y c <C 0.5) 
separated by narrow downflows with Y c w 0.49 (Fig. |l|). The shock reaches the Fe/Si interface 
at r = 1.4 x 10 s cm after ~ 100 ms. Shortly thereafter, at t w 120 ms temperatures right behind 
the shock have dropped below « 7 x 10 9 K, and 56 Ni starts to form in a narrow shell (Fig. |^). 
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Figure 2: Evolution of the 56 Ni mass fraction within the first 300 ms after the start of the calcu- 
lations. 56 Ni has only been synthesized in those regions which have an electron fraction Y e > 0.49 
and reached peak-temperatures above ~ 5 x 10 9 K. 



During the ongoing expansion and cooling, 56 Ni is also synthesized in the convective region. 
However, its synthesis proceeds exclusively in the narrow downflows which separate the rising 
bubbles and have a sufficiently high electron fraction Y c . This leads to a highly inhomogeneous 
nickel distribution (middle panel of Fig. ^) shortly before complete silicon burning freezes out at 
t ~ 250 ms. Moreover, convective motions are still present even after 56 Ni production ceases, and 
distort the nickel containing shell (upper panel of Fig. |2|). Only when convection stops around 
t ~ 400 ms, the flow pattern becomes frozen in and the entire post-shock region expands nearly 
uniformly. 

Due to the asphericity of the shock caused by the rising bubbles, "bent" shells containing the 
products of incomplete silicon burning as well as oxygen burning form outside the nickel-enriched 
region. The post-shock temperature drops below 2.8 x 10 9 K at t = 495 ms, when the shock is 
about to cross the Si/O interface. Thus, our model shows only moderate oxygen burning (due to 
a non- vanishing oxygen abundance in the silicon shell), and negli gible neon and carbon burning. 
This is caused by the specific structure of the progenitor model of Woosley ct al. (1985 ) and may 
change when different (especially more massive) progenitors are used. In total, 0.052 M Q of 56 Ni 
are produced, while 0.10 M© of material in the deleptonized bubbles are synthesized at conditions 
with Y c < 0.49 and end up as neutron-rich nuclei. The explosion energy of our 2D model saturates 
at 1.48 x 10 51 erg at t = 885 ms. This value is still to be corrected for the binding energy of the 
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Figure 3: Density distribution (log 10 P fe cm_3 ]) i n the exploding star, 100s after core bounce. 
The outer edge of the displayed domain (blue) coincides with a radius of 1.15 x 10 11 cm. The 
supernova shock (bright orange discontinuity near the outer edge) has crossed the He/H-interface 
and reached a radius of r = 10 11 cm. Rayleigh- Taylor mushrooms distribute the products of 
explosive nucleosynthesis throughout the He-core. The reverse shock (aspherical discontinuity) 
visible near the center was created after the main shock had passed the C+O/Hc-intcrface of the 
star. (The excessive outflows along the axis are specific to the use of polar coordinates.) 

outer envelope. 

4 Growth of Rayleigh- Taylor Instabilities 

During the next seconds the shock detaches from the formerly convective shell that carries the 
products of explosive nucleosynthesis, looses its asphericity and crosses the C+O/He-interface. 
Along with the slow-down which the shock experiences after passing this interface due to the 
varying density gradient, the entire post-shock material is also rapidly decelerated. Twenty seconds 
after core-bounce, the metal-containing shell has thus been compressed to a thin, dense layer which 
is bounded inwards by a reverse shock and contains two regions which show crossed density and 
pressure gradients. The first of these is located at the Ni+Si/O-interface while the second coincides 
with the C+O/He-interface. Thus, Rayleigh- Taylor instabilities at the Ni+Si/O and C+O/He- 
interfaces grow rapidly. 

At t = 100 s (Fig. ||) the instabilities are fully developed and have already interacted with each 
other. Nickel and silicon are dragged upward into the helium shell in rising mushrooms on angular 
scales from 1° to about 5°, whereas helium is mixed inward in bubbles. Oxygen and carbon, 
located in intermediate layers of the progenitor, are swept outward as well as inward in rising and 
sinking flows. At t = 300 s the densities between the dense mushrooms and the ambient medium 
differ by factors up to 5 while the fastest mushrooms have already propagated out to more than 
half the radius of the He core (Fig. [IJ) . 

We observe that the nickel is born with very high velocities, of the order of 15,000 km/s at 
t = 300 ms. These velocities drop significantly during the subsequent evolution, especially after 
the shock passes the C+O/He interface. At a time of t = 50 s most of the 56 Ni is expanding with 
3200 - 4500 km/s and maximum velocities Ujjf* are around 5800 km/s. During the following 50 
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Figure 4: Density distribution (log 10 p [gem -3 ]) 300 s after core bounce in the inner ~3x 10 11 cm 
of the star. The supernova shock (outermost blue discontinuity) is located inside the hydrogen 
envelope at r = 2.7 x 10 11 cm. A dense shell (visible as a bright red ring) has formed behind the 
shock. Its outer boundary coincides with the He/H interface while its inner boundary is in the 
process of steepening into a reverse shock. The fastest mushrooms have almost caught up with 
the inner boundary of the shell. (The excessive outflows along the axis are specific to the use of 
polar coordinates.) 

seconds, v™? x drops to ~ 5000 km/s. However, after t = 100 s, the clumps start to move essentially 
ballistically through the helium core and only a slight decrease of u^f* to ~ 4700 km/s occurs until 
300 s, when the bulk of 56 Ni has velocities below 3000 km/s. 

We have recently accomplished to follow the subsequent evolution of our model up to 16 000 s 
after core-bounce. The unsteady propagation speed of the supernova shock which has formerly 
induced the instability at the Ni+Si/O and C+O/He-interfaces also leads to the formation of a 
dense (Rayleigh- Taylor unstable) shell at the He/H interface. As before, the inner boundary of this 
shell steepens into a reverse shock (Fig. ||), while in the process of this the entire shell is rapidly 
slowed down. Our high-resolution simulations reveal a potentially severe problem for the mixing 
of heavy elements into the hydrogen envelope of Type II supernovae like SN 1987A. We find that 
the fast nickel containing clumps, after having penetrated through this reverse shock, dissipate a 
large fraction of their kinetic energy in bow shocks created by their supersonic motion through the 
shell medium. This leads to their deceleration to ~ 2000 km/s in our calculations. Contrary to all 
previous studies, which tried to accelerate the nickel by the instabilities, our computations show, 
that the main problem for obtaining high nickel velocities is how to avoid decelerating the clumps 
once they reach the He/H interface. The growth time scale of the instability at this interface is too 
large to allow for a fast shredding of the dense shell and the formation of "holes" through which 
the clumps could penetrate more easily. The high iron velocities in SN 1987 A might therefore 
hint towards a smoother density profile exterior to the He-core, which would suppress dense shell 
formation, or towards a global asymmetry of the explosion. 

During the computations we became aware of oscillations with angle in parts of the postshock 
flow (Figs. | and ^). These are caused by the "odd-even decoupling" phenomenon associated with 
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grid-aligned shocks (Quirk 1994; see also LeVeque 1998). As a consequence, the maximum nickel 
velocities, vS ax , obtained in our AMR calculations have probably been overestimated by ~ 25% 
because the growth of some of the mushrooms was influenced by the perturbations induced by 
this numerical defect. The main results of our study, however, are not affected. 



5 Conclusions 

We have studied the evolution of a core collapse supernova in a 15 M Q blue supergiant progenitor 
focusing on the first 300 s after core bounce. As a result of the interplay between neutrino-matter 
interactions, hydrodynamic instabilities and nucleosynthesis in the framework of the neutrino- 
driven explosion mechanism, the products of explosive silicon and oxygen burning are mixed 
throughout the inner half of the helium core of the progenitor star. Ballistically moving, metal- 
enriched clumps with velocities up to more than ~ 4000 km/s are observed. Rayleigh- Taylor 
instabilities at the C+O/He and Ni+Si/O-interfaces transport helium deep inward and sweep 
nickel, silicon, oxygen and carbon outward in rising mushrooms. 

Our simulations suggest that high- velocity metal-rich clumps are ejected during the explosion of 
Type lb (and Ic) supernovae. In case of Type II supernovae, however, the dense shell left behind by 
the shock passing the boundary between helium core and hydrogen envelope, causes a substantial 
deceleration of the clumps. Thus, we cannot confirm the results of Herant fc Benz (1992| ) that 



''premixing" of the 56 Ni within the first minutes of the explosion, and its later interaction with 
the instability at the He/H-interface, can explain the high iron velocities observed in SN 1987 A. 
Moreover, our calculations strongly indicate, that all computations of Rayleigh- Taylor mixing in 
Type II supernovae carried out so far (including the case of SN 1987 A) have been started from 
overly simplified initial conditions since they have neglected clump formation during the first 
minutes of the explosion. In future work, it will be most interesting to study how the discussed 
effects depend on different progenitor models and to explore the implications of the strong outward 
mixing of 56 Ni for the light curves and spectra of Type lb supernovae. 
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